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Abstract 

Pseudo-rapidity distributions of two high multiplicity events Ca-C and Si-AgBr 
observed by the JACEE are analysed by the wavelet transform. Wavelet spectra 
of those events are calculated and compared with the simulation calculations. The 
wavelet spectrum of Ca-C event somewhat resembles to that simulated with the 
uniform random numbers. That of Si-AgBr event is not reproduced by simulation 
calculations with Poisson random numbers, uniform random numbers, or a p- 
model. 

Introduction. In high energy ncleus-nucleus (AA) collisions, number density of sec- 
ondary particles in the rapidity space become very high, and studies of number density 
fluctuations in the rapidity space is expected to reveal new features of multiparticle 
production mechanisms. We have analysed pseudo-rapidity distributions in AA colli- 
sions including JACEE events, using Higuchi's method (a kind of length method), the 
auto-regressive model 0, and the fast Fourier transform ||. In order to classify high 
multiplicity events by those pseudo-rapidity distribution patterns, the fractal dimensions 
of the distributions (box dimensions) are estimated with both methods. 

In this paper, the pseudo-rapidity distributions of the two JACEE events, Ca-C and 
Si-AgBr are analysed by the wavelet transform 0,10, 0,0. Any function (data) can 
be expanded by self-similar wavepackets in this scheme. Therefore characteristics of 
local fluctuations can be extracted from the data. Wavelet spectra of the two events are 
calculated and are compared with simulation calculations. In the second paragraph, the 
wavelet transform is concisely introduced. In the third one, wavelet spectra of the two 
events are calculated. Those are compare with the simulation calculations in the next 
paragraph. The final paragraph is devoted to the concluding remarks. 
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Wavelet transform. Wavelets are constructed from dilation and translation of a scal- 
ing function. The scaling function <f)(x) is constructed by an iteration equation, 



N-l 

(x) = E c fc &_i(2a; - fc) 1=1,2,-.., (1) 

fc=0 



from a primary scaling function (f>o(x), where N is a even number, and Ck (k — 0, 1, • • • , N— 
1) are constants. The iteration is continued until 4>i{x) becomes indistinguishable from 
4>i-i(x), and (j)(x) is defined by 4>i{x). The primary scaling function 4>q(x) is taken as 



1 for < x < 1, 

otherwise. 



(f) (x) = 

The mother wavelet W(x) is given by 

N-l 

W(x)=Y,(- l ) k Ck<P(2x + k-N+l), (2) 

k=0 

and the j-th level wavelet (j — 0, 1, • • •) is defined as 

V^(z) = 2^(2^ - fc), fc = 0, 1, • • • , V - 1. (3) 

Coefficients c& (k — 0, 1, • • • , N — 1) should be determined so that the wavelets and the 
scaling function satisfy the following orthogonal conditions, 

J 4>{x) (ft(x)dx = 1, 

J <f>(x) ^j, k (x)dx = 0, (4) 

If N=4, the coefficients are given by 

l + v^ 3 + V% 

CO = — Cx = — 

(5) 

3 - 1 — >/3 

Arbitrary function /(#) defined in the region < x < 1 can be expanded as 

/ = /* + E/ w » 

/* = a o 0(x), (6) 

2^-1 

/0) = E Oij, k ^j,k{ X )- 
k=0 
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By the use of Eqs.(|) and (|5|), We have a relation, 

J f{xfdx = a 2 + Y, E ji 



(7) 



3=0 



where Ej denotes the j-th level wavelet spectrum 



E 



23-1 

E 



/ / a j,ki 
k=0 



J = 0,1, 



Equation (|7|) corresponds to the Parseval's equation in the Fourier transform. 

Hereafter we consider the case that the function f(x) defined in the region < x < 1 
is given by descrete value, 



f i = f{x 1 + A-(i-l)), 



1,2,- --,2 r . 



We assume that f(x) is periodic outside the defined region. As the data points are 2 r 
f(x) is expanded by finite terms. Summation on the right hand side (RHS) of the Eq.([5 
runs from j = to j = r — 1. 

By the use of column matrices, 



/ 



/a 
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Eq.( 



is expressed as 
/ - 

tj(r-l) . 



V a j,2i-l J 



3=0,1, 



25 a LiL 2 • • -L r , 

2 2 Aj_xHjLj + i ■ ■ ■ L r , 

2 5 ' A r _\H r , 



i = 1,2, 



(9) 



where Lj and ifj are 2 J 1 x 2 J matrices. The z— th row of Lj has 2(z — l)'s zeroes from 
the first to the 2(i — 1) -th elements, and is expressed as, 



(° 







2(i-l)+l *2(i-l)+2 



v/2 



c ci 



2(i-l)+7V 







The z— th row of iJj also have 2(i — l)'s zeroes from the first to 2(i — l)-th elements and 
is written as 



(° 



' ^2(i-l)+l "-2(i-l)+2 
1 

71 







Cat_i -Cat- 



h-2(i-l)+N 



c 
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If 2(i — 1) + k > 2 J in Lj ( or Hj ), element Z 2 (i-i)+fc ( or /i 2 (i-i)+fc ) is added to the 
(2(z — 1) + A; — 2 J ) -th element in each row. 

From Eq.(|5|), matrices Lj and Hj satisfy the following conditions, 

L j L j = I , Z/j £^ = , 

(10) 

£j *_Lj = 0, Hj ^ Hj = /, 

where t Lj denotes the transpose of matrix Lj, and I is the 2 J ~ 1 — th order unit matrix. 
Then, the wavelet coefficients are written as 

r 

ao = 2 2 LiL 2 ■ ■ • L r f, 
Aj-! = 2- r *H j L j+1 ---L r f : j = 1,2,-.., r-l, (11) 

As the summation on the RHS of Eq.(^) runs from j = to j — r — 1 in this case, the 
wavelet spectra Ej from j = to j = r — 1 are obtained. 

Wavelet spectra of the data. Pseudo-rapidity rj distributions of Ca-C and Si-AgBr 
events are shown in Fig. la and b, respectively. From each distribution, we subtract the 
background distribution 0, 

f ( V ) =A[(1- e- y -")(l - e-^)] B , —Y <r]<Y, (12) 
4 = 81, £ = 5.4, y = 7.0 for Ca-C, 
A= 184, £ = 8.1, Y = 5.5 for Si-AgBr. 



Pseudo-rapidty distributions of Ca-C and Si-AgBr events are given with the bin size 
Arj = 0.1, and the number n of the data points in each event is in the range 2 6 < n < 2 7 . 
Those used in the analysis should be 2 r . Therefore 64 points (r = 6) are used for both 
events; —3.0 < i] < 3.4 for Ca-C, and —2.7 < rj < 3.7 for Si-AgBr. Fig 2a and b show 
the distributions of the two events, where the backgrounds are subtracted. The standard 
deviation a of the distribution shown in Fig.2a is 2.71, and that of fig. 2b is 4.03. 

In Fig. 3a and 3b, the wavelet spectra Ej (J = 0, 1, • • ■ , 5) of Ca-C and Si-AgBr are 
shown, respectively. The black circles show the wavelet spectra of the distributions 
where backgrounds are subtracted. The blank circles show that of the raw distributions 
( without subtraction of the backgrounds). The results are different at the 0-th and the 
first level, but those from j =2 to j=5 are almost the same. The wavelet spectrum of Ca-C 
event increases linearly in semi-logarithmic scale from j=2 to j=5. That of Si-AgBr is 
approximately flat in the same range. 

Comparison with simulation calculations. The wavelet spectra of the two events are 
compared with the following simulation calculations; 



5 



(i) Poisson random numbers with a mean \x. The standard deviation a of the Poisson 
distribution with the mean fi is given by a = ^JJi. 

(ii) uniform random numbers from to a. The mean \x is given by, 

\i — - \ xdx =—. 
a Jo 2 

The variance a 2 is given by 

a 2 = — [ (x — n) 2 dx = — . 
a Jo K 12 

Then the parameter a is expressed by the standard deviation as 

a = 2Vocr. 



(iii) a p-model|10[ with a fraction p a (p& = 1 — p a ) and an initial 'energy' (or particle 
number) Eq. 

In the p-model, the initial energy E in a rapidity interval Ay is devided into p a E 
and pb E Q at the first step, as the interval Ay splits into two sub-intervals with 
equal width. p a E is regarded as the energy in one interval, and pb E is that in 
the other interval. At the subsequent step, each energy is devided into two in the 
same way, as each sub-interval splits into two. After n steps, the initial energy E 
is devided into 2 n terms; energy p ] a Pb~ J E (j = 0,1, • • • ,n) appears n Cj times. 
Those are energies in 2 n ordered sub-intervals of Ay with the same width. The 
mean /j, and the variance a 2 are given respectively by 

1 n E 

Z 3=0 Z 

^ = ^ t (ftrftf-O'-Q -/•"= (f)' (2- (rf + £)"-!). 

In the simulation calculations, we iterate r + 2 = 8 steps, and randomize the 
ordering of the energies. Then, we add every four energies, and have 2 6 ones, 
which are compared with the data. 

Each simulation calculations are done 1000 events. The mean value is not sensitive 
to the wavelet spectra, and parameters are adjusted to reproduce the standard deviation 
of each distribution shown in Fig. 2a or b. 

We calculate the wavelet spectra Ej imu (j = 0, l,---,5) of simulated events and 
count the number of events within a value of z, 

_ 5 

£(log 2 £,-log 2 £f™) 2 . (13) 

J=2 
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In order to generate the random numbers, we use a congruence method[[Tl]]. At first, 
we choose an integer r , using the equation 



■ a Tj^i + c 
a = 1229, 



mod m 



(14) 



c = 351750 



m = 1664501 



we generate the random numbers rj (j = 1,2, •••) subsequently. From those random 
numbers we get the uniform random number 



The results for Ca-C and Si-AgBr events are shown in Table la and lb, respectively. 
Table la shows that the wavelet spectrum of simulated events for Ca-C events have more 
than 15 % out of 1000 simulated events within z < 0.8 in all the three cases. Those with 
uniform random numbers have 20 % and is higher than the other two cases. 

Contrary to the Ca-C event, simulated events for Si-AgBr event have scarecely low 
rate within z < 0.8 in all the three cases (See table lb). The wavelet spectrum of Si-AgBr 
event cannot be reproduced by any of these simulation calculations. 

Examples of the simulation calculations with (i) Poisson random numbers are shown 
in Fig. 4a and b, with (ii) uniform random numbers in Fig.5a and b, and with (iii) the 
p-model in Fig 6a and b. 

Concluding Remarks. Two high multiplicity events Ca-C and Si-AgBr are analysed 
by the wavelet transform. Wavelet spectra of the two events are calculated and com- 
pared with the simulation calculations. Wavelet spectrum of the Ca-C event somewhat 
resembles to the simulation calculations with uniform random numbers, but that of the 
Si-AgBr event can not be reproduced by any of the three simulation calculations. Further 
investigation will be reported elsewhere that there are any simulation calculations which 
reproduce n distributins of high multiplicity events. Our analysis show that observed 
high multplicity events can be classified by the wavelet spectra, and what statistical 
regularity (or irregularity) would be dominant event by event. 
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s i = — 
m 



< Si < 1. 
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Table Caption 



Table 1 Number of simulated events which satisfy the condition z < 0.4, z < 0.6, or z < 0.8 
for a) Ca-C and b) Si-AgBr. Each simulation calculations are done 1000 events. 



Figure Captions 

Fig.l Pseudo-rapidity distributions of a) Ca-C and b) Si-AgBr events. 

Fig. 2 Distributions used for wavelet analysis. Those are obtained from rj distributions 
where the backgrounds are subtracted; a) Ca-C, and b) Si-AgBr. 

Fig. 3 Wavelet spectra of a) Ca-C and b) Si-AgBr events. Black circles are calculated 
from the rj distributions where the backgrounds are subtracted. Blanck circles are 
from rj distributions ( where backgrounds are not subtracted). 

Fig. 4 Wavelet spectra of simulated events with Poisson random numbers for a) Ca-C 
with r = 44255 (z=0.14) and b) Si-AgBr with r = 115383 (z=0.50). 

Fig. 5 Wavelet spectra of simulated events with uniform random numbers for a) Ca-C 
r = 6634423 (z=0.25) and b) Si-AgBr with r = 115383 (z=0.55). 

Fig. 6 Wavelet spectra of simulated events with the p-model for a) Ca-C with r = 
1567321 (z=0.27) and b) Si-AgBr with r = 707048 (z=0.58). 
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This figure "figl-l.png" is available in "png" format from: 
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